Semiparametric Regression in Size-Biased Sampling

Summary .  Size-biased sampling arises when a positive-valued outcome variable is sampled with selection probability proportional to its size. In this article, we propose a semiparametric linear regression model to analyze size-biased outcomes. In our proposed model, the regression parameters of covariates are of major interest, while the distribution of random errors is unspecified. Under the proposed model, we discover that regression parameters are invariant regardless of size-biased sampling. Following this invariance property, we develop a simple estimation procedure for inferences. Our proposed methods are evaluated in simulation studies and applied to two real data analyses.     

Spatial and temporal factors associated with nest survival of Gray Flycatchers in managed ponderosa pine forests

Gray Flycatchers (Empidonax wrightii) breed in a variety of habitats in the arid and semi‐arid regions of the western United States, but little is known about their breeding biology, especially in the northern portion of their range where they nest in ponderosa pine (Pinus ponderosa) forests. From May to July 2014 and 2015, we conducted surveys for singing male Gray Flycatchers along the eastern slope of the Cascade Range in Washington, U.S.A, monitored flycatcher nests, and quantified nest‐site vegetation. We used a logistic‐exposure model fit within a Bayesian framework to model the daily survival probability of flycatcher nests. During the 2 yr of our study, we monitored 141 nests, with 93% in ponderosa pines. Mean clutch size was 3.6 eggs and the mean number of young fledged per nest was 3.2. Predation accounted for 90% of failed nests. We found a positive association between daily nest survival and both nest height and distance of nest substrates from the nearest tree. Flycatchers that locate their nests higher above the ground and further from adjacent trees may be choosing the safest alternative because higher nests may be less exposed to terrestrial predators and nests in trees that are farther from other trees may be less exposed to arboreal predators such as jays (Corvidae) that may forage in patches with connected canopies. Nests in trees farther from other trees may also allow earlier detection of approaching predators and thus aid in nest defense.     

Simulation of Y-chromosomal haplotype data

The non-recombining nature of the Y-chromosome determines the non-independence of alleles between loci. The evolution of short tandem repeat (STR) loci in the Y-chromosome is the result of different factors such as differential mutation rates, mutation modes, gene conversion, selection and demographic processes. The degree of correlation between loci is dependent on the magnitude of these processes. The simulation of data is a routine tool used for testing hypotheses in population and evolutionary studies. The most basic parameters hitherto used in lineage haplotype simulations are the allele frequency distributions and mutation rates, assuming either full independence or linkage between loci. In this study we introduce use of the Spearman correlation coefficient to estimate the degree of dependence between non-recombining loci. Then, both the interdependence between loci and the allele frequency distributions at multi-allelic loci are incorporated in an algorithm for simulating haplotypes. We illustrate the method using published and unpublished Y-chromosome STR data.     

Can transposable element copy number distribution parameters be estimated from natural populations of Drosophila melanogaster?

The copy frequency distribution of a transposable element family in a Drosophila melanogaster natural population is generally characterised by the values of the Charlesworths' model parameters α and β (Charlesworth & Charlesworth, 1983). The estimation of these parameters is made using the observed distribution of the occupied sites in a population sample. Several results have been interpreted as due either to the influence of stochastic factors or to deterministic factors (transposition, excision, selection…). The accuracy of this method was tested by estimations performed on samples from simulated populations. The results show that with the sample size usually used for natural population studies, the confidence intervals are too large to reasonably deduce either the element copy number distribution or the values of transposition and excision rate and selective coefficients.     

Joint multiple quantitative trait loci mapping for allometries of body compositions and metabolic traits to body weights in broiler

In order to map quantitative trait loci (QTLs) for allometries of body compositions and metabolic traits in chicken, we phenotypically characterize the allometric growths of multiple body components and metabolic traits relative to BWs using joint allometric scaling models and then establish random regression models (RRMs) to fit genetic effects of markers and minor polygenes derived from the pedigree on the allometric scalings. Prior to statistically inferring the QTLs for the allometric scalings by solving the RRMs, the LASSO technique is adopted to rapidly shrink most of marker genetic effects to zero. Computer simulation analysis confirms the reliability and adaptability of the so-called LASSO-RRM mapping method. In the F₂ population constructed by multiple families, we formulate two joint allometric scaling models of body compositions and metabolic traits, in which six of nine body compositions are tested as significant, while six of eight metabolic traits are as significant. For body compositions, a total of 14 QTLs, of which 9 dominant, were detected to be associated with the allometric scalings of drumstick, fat, heart, shank, liver and spleen to BWs; while for metabolic traits, a total of 19 QTLs also including 9 dominant be responsible for the allometries of T4, IGFI, IGFII, GLC, INS, IGR to BWs. The detectable QTLs or highly linked markers can be used to regulate relative growths of the body components and metabolic traits to BWs in marker-assisted breeding of chickens.     

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W1 and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the O_b atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.     

A Molecular Dynamics Investigation of Chain Conformational Changes In Compressed Bilayers

Abstract

The conformations of the chains constituting the hydrophilic component of alkyl monolayers and bilayers are investigated by performing molecular dynamics atomistic simulations on these systems at different temperatures. Several monitoring techniques are used to reveal the chain conformations, including atom pair radial distribution functions, evolutions of the torsional angles over thousands of timesteps, frequency distributions of the torsionl angles and ‘snapshot’ plots of the atomic configurations. These methods consistently testify to the stability of the trans (fully extended) character of the strain-free alkyl chains up to room temperature. The chains retain much of this conformation even when the layers are compressed by the application of pressure, to which the chains respond by ‘folding’ at the ends attaching them to the substrate planes while maintaining directions which are mainly normal to these planes. A non-zero gap between the layers is also maintained. A pressure of about 50 kbar abruptly causes all motion in the chains to cease, resulting in a highly ordered lattice structure.